The absolute position of a resonance peak 
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It is common practice in scattering theory to correlate between the position of a resonance peak 
in the cross section and the real part of a complex energy of a pole of the scattering amplitude. 
In this work we show that the resonance peak position appears at the absolute value of the pole's 
complex energy rather than its real part. We further demonstrate that a local theory of resonances 
can still be used even in cases previously thought impossible. 
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There are few phenomena in physics which are as strik- 
ing and ubiquitous as resonance peaks in the cross sec- 
tion. As such, it is of great importance to understand 
their origin and to be able to predict the resulting reso- 
nance profile in the cross section. Over the years many 
methods have been proposed and successfully used to an- 
alyze and explain resonance peaks appearing in the cross 
section. Among which we mention the use of the poles of 
the scattering matrix(S- matrix) [ij, the Fano configura- 
tion interaction y, and the Feshbach resonance theory 

i- 

The main purpose of such theories is to reconstruct the 
resonance structure appearing in the cross section. This 
should be done with but a few parameters which can be 
used to characterize the scattering system at hand. It is 
in the way these parameters are chosen and calculated 
that the methods mentioned differ from each other. One 
such difference is the use of non-local versus local pa- 
rameters, i.e., whether the system parameters used to 
describe the resonance profile are energy dependent or 
not. 

The most famous application of a local resonance the- 
ory is the Breit-Wigner(BW) resonance profile [4| which 
connects an appearing resonance peak in the cross sec- 
tion with a pole of the S-matrix in the complex energy 
plane. The resulting resonance profile reads: 



above equation for an arbitrary scattering potential. 
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where, e is known as the resonance position and is given 
by the real part of the pole in the energy plane, and F is 
width of the Lorentzian at half maximum which is twice 
the imaginary part of the resonance pole in the energy 
plane. Although the above resonance profile has a very 
impressive track record it fails whenever the pole is not 
isolated (from the effect of other poles), it is near the 
threshold, or it is far from the real axis. One possible 
resolution to this predicament is to move into a non- 
local description and define and energy dependent width 
and position such that the resonance profile now reads. 
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Although clearly able to 



■^ 1 


i<* 


^ ^ ^ «k 


sion Probabilit> 

o o 
'(J) bo 


r 


En 


\i 


En\ 


2.06- 2. 15z 


3 


2.97 


■PO.4 




10.72 - 5.25z 


12 


11.94 


E 


/ 


21.41 - 8.20z 


23 


22.93 


20.2 

H 

n 


/ 





10 15 

Energy 



20 



25 



reproduce any resonance structure it is neither straight 
forward nor simple to calculate the non-local terms in the 



FIG. 1. (color online) The transmission probability above 
a square well potential of depth Vo = — 13 and length L = 
7r/\/2. Full(green) and dashed(red) arrows mark the real part 
and absolute value of the resonance poles energy respectively. 
The table in the inset shows the values of the first three en- 
ergies, En, of the resonance poles, the analytical position, 
A„, of the transmission unities, and the absolute value of the 
resonance energies. 

Recently, the use of the poles of the scattering matrix 
was reinvigorated by the novel method presented by Tol- 
stikhin et al. [7[ . The possibility to solve the Schrodinger 
equation with Siegert boundary conditions [8| to obtain 
all the poles of the scattering matrix opened the door 
to describe the entire cross section using only the poles. 
Tracing the origin of the structures in the cross section to 
specific poles of the matrix has always been invaluable in 
order to understand the scattering mechanisms respon- 
sible for the resulting cross section. It would, therefore, 
be extremely beneficial if one could assign to every reso- 
nance in the cross section a specific pole. It is the failure 
of this possibility that necessitates the use of non-local 
theories and as we shall demonstrate there is still much 
to be done using only local resonance profiles. 



In order to demystify the above discussion let us start 
with the simplest example where one cannot connect be- 
tween the poles of the S-matrix and the resonance ap- 
pearing in the cross section using the assumptions of the 
BW profile. Consider a one dimensional square well. It 
is well known, see for example [5|, that the transmis- 
sion spectrum displays an oscillating behavior and the 
transmission probability reaches unit value only at dis- 
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Crete energies satisfying: A„ = — Vb + 2mLA ' ^^'-^ ^^^^ 
T{Xn) — 1. In the above equation Vq and L are the 
well's depth and length respectively, m is the mass of 
the projectile, and n is a positive integer. Fig. [1] de- 
picts the transmission for a specific choice of the wells 
depth and length. The poles of the scattering ampli- 
tude can be calculated by solving the Schrodinger equa- 
tion with outgoing boundary conditions [6|. The posi- 
tions of the first three resonance poles, i.e., real part of 
the complex energy of the poles appearing in the fourth 
quadrant, are depicted in Fig. [T]using full(green) arrows. 
As can be readily observed one cannot directly assign a 
transmission unity with the real part of the resonance 
pole's energy. At first this may seem obvious since the 
resonance poles clearly overlap no direct connection be- 
tween a single pole and the observed resonance can be 
expected. But, as we shall show in this letter there is 
a one to one correspondence between each transmission 
unity and a single resonance pole. To relieve the sus- 
pense we give here the final answer that the resonance 
positions in the cross section are related to the absolute 
value of the pole's complex energy rather than to its real 
part. This is demonstrated in Fig. [T]by the dashed(red) 
arrows indicating the absolute values of the complex res- 
onance poles' energies. 

Let us consider a scattering process from a sin- 
gle one dimensional channel. As mentioned above, 
the transmission amplitude can be expressed using 
all the poles of the scattering matrix [7[. These 
are found as solutions of the Schrodinger equation : 

[-^dl + V{x)]'i'n{x) = ^"ifnix) along with Siegert 
boundary conditions imposed at x — ±L which satisfy 
[i|: ((9^*„(±L)=F«fcn*«(±-^)) = 0. Here we shah use 
atomic units for which h = m — I. We choose to nor- 
malize the Siegert states such that: 2ifc„ J_^ '^'^{x)dx — 
[^liL) + *.^j(-L)] =1. For all but piece- wise poten- 
tials the solutions must converge with respect to L. One 
can categorize the Siegert solutions according to where 
they fall in the complex k-plane: Bound(anti-bound) 
solutions on the positive(negative) imaginary axis and 
resonance(anti-resonance) solutions on the fourth(third) 
quadrant. Using the relation En — kfj^, the solutions' 
energy can be found. 

Using the Siegert states solutions, the transmission 
amplitude for a one dimensional potential can be writ- 
ten as 0: t = 2ke-^^>^^ En=i '^"^ti-"fc"'^^ ^^lere the 
sum runs over all poles, i.e., bound, anti-bound, virtual, 



and resonance poles, and k = \J2E. The sum above can 
be somewhat simplified by transforming it into a prod- 
uct. Such a product formula was first given by Tolstikhin 
et al. in [lO| for symmetric potentials and by Ostro- 
vsky and Elander in ll| for asymmetric potentials. For 
symmetric potentials, the product formula reads: t = 

lln=l 
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k-kZ, 



where N^{NJ) 



correspond to symmetric(anti-symmetric) Siegert states. 
Since the poles fc„ are either purely imaginary or come 
in pairs fc„ and — fc*, the above product terms can each 
be written as a pure phase, the well known Blaschke 
factor |12|. Considering first the pairs of resonance 
and virtual states one can reformulate the product 



of each pair as: 



fc+fc^ fc-(fc„)'' 
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where after 



some trigonometric manipulations, one finds that S^ 
( k (r±/2) • 



arctan I 



the resonance energy: 



with the common definition of 



I^n — 9 



£„ — i-^. For 
the bound and anti-bound poles each term in the prod- 
ucts read: (k + kn)/{k — k^) = e^*"^" , where for these 
poles 6n = arctan (Im[fc„]/fc). Defining A± = J^ri^i^n 
where the sum runs over all bound poles, anti-bound 
poles, and resonance-virtual pole pairs, the transmis- 
sion amplitude can be brought to the following form: 
t = ie-2i'=ie*(^++^-)sin(A+ - A_), and the transmis- 
sion probability reads, T = sin^(A+ — A_). 

1 



^0-8 

O 

qI 0.6 

c 
o 

CO ^-^ 

E 

c 0.2 
CD 



: ] '. \\ A 2:^ I 

; J \ »\'-3 -2-1 12 3 

/ V/^^---, 



0.1 



0.2 0.3 

E [a.u.] 



0.4 



0.5 



FIG. 2. (color online) The numerically exact transmission 
probability (solid(black) line) through a double barrier po- 
tential, see the text, using the following parameters: /3 = 5/2, 
7 = 0.8, and a — 0.5. We show only the first resonance peak. 
On top of the exact resonance peak we draw the BW reso- 
nance profile (dotted(green) line) (see Eq. [l| and the novel 
resonance profile given here (dashed(red)line) (see Eq. O) 

Using the above formula for the transmission proba- 
bility, one can examine the contribution of a single reso- 
nance pole to the cross section. Without loss of general- 
ity, let us assume that the corresponding Siegert eigen- 
state is even. Then neglecting the contribution of all 



but a single resonance-virtual pole pair, the transmission 
probability reads, 



T: 
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where the resonance pole is located at kr, 

kL. 
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trigonometric identity: sin(arctan(a:)) = 
sight, Eq. [21 is very reminiscent of the BW resonance 
profile, see Eq. [TJ There are, however, several crucial 
differences. First and foremost, the resonance peak has 
a maximum at the absolute value of the complex reso- 
nance pole energy. We have already demonstrated above, 
as shown in Fig. [TJ that this allows for the connection 
of each transmission unity with a single resonance pole 
even in the case of wide overlapping resonances. The 
resonance profile in Eq. [5] also differs from the BW pro- 
file in that it is not a Lorentzian and is asymmetric with 
respect to its maximum, rising sharply to the left of the 
maximum and falling more slowly to the right of the max- 
imum. Even though extension of the BW profile into an 
asymmetric one are known [l|, the current profile for- 
mula characterizes the asymmetry with only the reso- 
nance complex energy, i.e., with the same information 
needed to construct the BW profile, and does not require 
the calculation of additional parameters. 

Wide overlapping resonances are not the only case 
where the local description of resonance is abandoned 
for a non-local one. For resonances near the threshold, 
a similar belief that one cannot correlate between a res- 
onance peak and a single pole exists. However, the reso- 
nance profile developed here shows that even very close 
to the threshold the resonance peak is due to a single 
resonance pole, and one can reconstruct the resonance 
peak in the cross section using only the resonance pole 
energy. Fig. [21 depicts the transmission probability near 
an isolated resonance in a double barrier potential of the 
form: V{x) — (^(3x^ — 7) e^°'^ , see the inset. Clearly, the 
profile in Eq. [2] reproduces the peak in the transmission 
probability very well using the same information used in 
the BW resonance profile which as seen in the figure fails 
to reproduce the resonance peak quantitatively. Thus 
one does not need to resort to non-local theories to de- 
scribe resonance peaks near the threshold even when, as 
is in the case portrayed, the imaginary part and real part 
of the resonance energy are comparable in magnitude. 

We turn now to discuss yet another important case 
where the local approximation has be erroneously dis- 
carded. In his seminal paper Nussenzveig [13| studied 
the effect of varying the potential on the poles of the 
scattering matrix. One of the most studied aspects has 
been the behavior of the poles as the potential is varied 
to support an additional bound state. It is well estab- 
lished that for a single s-wave like channel, looking in 
the k-plane, the resonance poles located in the fourth 
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FIG. 3. (color online) The first resonance peak in the trans- 
mission probability for the double barrier potential, see the 
text, for different values of 7: 0.875(a), 0.885(b), and 0.89(c). 
The numerically exact transmission (solid(black) line) is com- 
pared for each value of 7 with the novel resonance profiles 
presented here in Eqs. [2] and [3] (dashed(red) line). The in- 
set portrays the poles of the scattering matrix closest to the 
imaginary k axis for each of the cases (a), (b), and (c). For 
convenience we also mark the 45 degree bisector indicating 
in the third and fourth quadrants. Poles in the fourth(third) 
quadrant to the left(right) of the bisector have a negative po- 
sition. 



quadrant and the virtual state poles located in the third 
quadrant move toward the negative imaginary k-axis as 
the potential well is, for example, deepened. Eventually, 
two poles approaching from either side of the imaginary 
axis coalesce, thereafter, forming two anti-bound state 
poles which lie on the negative imaginary axis. Follow- 
ing the creation of two anti-bound poles, one of these 
poles moves up toward the origin and eventually leads 
to the formation of a new bound state, while the other 
moves down on the imaginary axis. 

The connection between the motion of the resonance 
pole described above and the corresponding motion of 
the resonance peak in the cross section remains some- 
what elusive even today. The main problem that arises 
in the analysis is that at some point, while approaching 
the negative imaginary axis, the resonance pole crosses 
the 45 degree bisector of the fourth quadrant of the k- 
plane making the resonance position negative. The reso- 
nance peak in the cross section, however, is still observ- 
able and is clearly situated above the threshold. Nussen- 
zveig deemed these resonance poles with negative posi- 
tion as the only "meaningless" [13| poles. As we shall 
show they are just as meaningful as any other pole of the 
scattering matrix. 

Consider the effect of deepening the well between the 
two barriers, see the potential depicted in Fig. [21 on the 
transmission spectrum. At some critical well depth, a 



new bound state appears. We wish to study the behav- 
ior of the resonance peak in the transmission spectrum 
as the potential depth is increased towards this critical 
value. Previous studies noted that as the well deepens 
the resonance peak approaches the origin. Figure [3] por- 
trays the first resonance peak in the transmission spectra 
for three different depths' of the potential well. As can 
be seen in the figure, the resonance peak moves toward 
the origin as we approach the critical value for which a 
new bound state is formed. Aside from the obvious shift 
in position, the different resonance peaks look almost the 
same. This is rather surprising since as we shall explain in 
the following, the potential responsible for the left most 
peak (marked (c)) doesn't support a resonance pole in 
the vicinity of the position of the peak. 

Let us examine first the two peaks in figure [3] marked 
(a) and (b). The potential responsible for peak (a) sup- 
ports a resonance pole above the threshold similar to that 
depicted in figure [5J As can be seen the resonance profile 
given in Eq. [2] (dashed(red) line) reproduces the numeri- 
cally exact resonance peak in the spectrum also depicted 
(solid(black) line). For the peak marked (b), however, the 
potential is such that the resonance pole moved passed 
the 45-degree bisector of the fourth quadrant of the k- 
plane, see inset in figure [3] Therefore, the position of 
the resonances pole is now negative. As the figure shows, 
the profile in Eq. [2] can be used very successfully to 
reproduce the resonance peak using only the resonance 
pole with the negative position. This puts the resonance 
poles with negative positions on the same footing as other 
resonance poles of the scattering matrix. The resonance 
pole's trajectory in the complex k-plane brings about yet 
another difficulty. As the pole approaches the negative 
imaginary axis, the resonance width defined as F = 2kikr 
goes to zero since fc^ — >■ 0. One might therefore conclude 
that such narrow resonances will be very hard to ob- 
serve. Such a misconception stems from the expectation 
of a B W profile whose full width half maximum (FWHM) 
is F. The profile presented in Eq. [21 however, is asym- 

For over 



metric and has a FWHM of FW 1 + 2 

the threshold, narrow resonances this is approximately 
F. However, for negative position resonances close to the 
negative imaginary axis one can no longer neglect the 
change in the FWHM. For the peak marked (b) in fig- 
ure [31 the resonance pole width is of the order of 10^'^ 
whereas the peak appearing in the figure has a FWHM 
of the order of 10^^. We reemphasize that the resonance 
peak, position and FWHM, are characterized using only 
the resonance pole complex energy, i.e., the same infor- 
mation used in the BW profile. 

Further increase in the depth of the well, past the 
branch point where the two poles coalesced, produces two 
anti-bound states. Before proceeding, we point out that 
even though after the coalescence of the resonance and 
virtual poles, i.e., the resonance pole no longer exists, the 



peak is still apparent in the transmission spectrum. As 
figure [3] clearly shows, the transmission spectrum seems 
somewhat oblivious to the identity crisis the poles are ex- 
periencing. The change from resonance and virtual poles 
to two anti-bound states left no impression on the result- 
ing spectrum. Since the colliding resonance and virtual 
states had the same spatial symmetry, the resulting two 
anti-bound states also share the same symmetry. We can 
therefore examine the transmission spectrum due to two 
anti-bound states of the same symmetry. Neglecting all 
contribution but those of the two anti-bound (AB) poles 



located at k\ 



AB 



-ik^ and fc^'^ = — Jfc^, the transmis- 



sion probability reads 
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The remarkable agreement between the above formula 
and the numerically exact transmission can be witnessed 
in Fig. [21 

In summary, we have demonstrated that a local theory 
of resonances can still be used even in cases previously 
thought impossible. Our analysis allows one to unequiv- 
ocally connect the resonance peak in the transmission to 
specific poles of the scattering matrix throughout vari- 
ous changes of the scattering potential. Even for wide 
overlapping resonances one can connect each transmis- 
sion unity with a single resonance pole. 
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